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Abstract 

We solve numerically the Schwinger-Dyson (SD hereafter) ghost equa- 
tion in the Landau gauge for a given gluon propagator finite at fc = 
(cigiuon = 1) and with the usual assumption of constancy of the ghost- 
gluon vertex ; we show that there exist two possible types of ghost dressing 
function solutions, as we have previously inferred from analytical consid- 
erations : one singular at zero momentum, satisfying the familiar relation 
ctgiuon+'2aghost = betwceu the infrared exponents of the gluon and ghost 
dressing functions(in short, respectively qg and ap ^), and having there- 
fore Oghost ~ —1/2, and another which is finite at the origin {aghost ~ 0) , 
which violates the relation. It is most important that the type of solution 
which is realized depends on the value of the coupling constant. There are 
regular ones for any coupling below some value, while there is only one 
singular solution, obtained only at a critical value of the coupling. For all 
momenta k < 1.5 GeV where they can be trusted, our lattice data exclude 
neatly the singular one, and agree very well with the regular solution we 
obtain at a coupling constant compatible with the bare lattice value. 



1 Introduction 

Since the first attempts, an impressive progress has been made in understanding 
the solutions to the Schwinger-Dyson equations for the QCD propagators and 

^Laboratoire associe au CNRS, UMR 8627 

^Let us recall that we denote the gluon by a G, and then the ghost by a F, for "fantomes" 
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their behaviour at smaU momenta. In particular, an important step has been 
accomplished by putting forward the essential contribution of internal ghost 
loops in the gluon propagator equation, previously neglected ; it has been shown 
that it may completely change the previously expected behavior of the gluon 
propagator from much more singular than the free one ( something like 1/fc'*- as 
was believed for a long time), to being much less singular pj ^. All the following 
considerations are assuming the choice of the Landau (i.e. Lorentz) gauge. 

-The consensus before 2005 

For some years, a consensus seemed to be obtained around a statement of 1) 
an infrared fixed point of the gluon-ghost coupling 2) a singular ghost dressing 
function (see below for more explanation). This consensus was very strong and 
unopposed, since several other approaches were apparently converging to the 
same conclusions. Such authoritative people as S. Brodsky have also appealed 
to it to support their considerations about AdS/CFT ( see [2]**). Another part 
of the consensus, deduced from a solution to coupled SD equations, was the 
statement that aa > 1, i.e. the gluon propagator should necessarily vanish; yet 
this was contested by a thorough calculation of Bloch [S] . 

To be more specific about this consensus, a usual assumption for the in- 
frared behaviour of gluon and ghost dressing functions is that they should be 
power behaved, i.e. for the gluon G{k) and the ghost F(k) dressing functions 
respectively : 

G(fc) ^ {k^r- (1) 

F{k) - (fc2)"^ (2) 

In a number of studies [1] |3], it has been stated that for a suitably simple 
assumption concerning the ghost-gluon and gluon vertices, the Dyson-Schwinger 
(SD hereafter) coupled equations for G{k) and F{k) imply 

aa + 2aF^0 (3) 

This statement is the starting point for the popular claim of an infrared fixed 
point for the QCD renormalised coupling constant. In fact, admitting the va- 
lidity of eq. ([3]), the IR fixed point would be present in the coupling constant 
defined by the ghost-gluon 3 points Green function in a MOM scheme. Let us 
recall the renormalisation conventions, with bare quantities denoted by a "B" 
subindex. In general : 

Gsik^) = Zs Gnik''), Fsik'') - Z3 Fni^), gs = Zg gn, 

Tr - Ji Fs, Zg = {Zl'^ Z^)-^ (4) 

In the MOM schemes : 

^3 = Gb(m') (5) 

Z, = Fb{ii^) (6) 

Gfl(fc2,M) = GB(fc')/GB(M') (7) 
Fn{k\^l)^FB{e)/FB{^l'') (8) 

^For this particular case, see especially the section 3 of the quoted paper 
*See especially the section 3 of the paper by Brodsky giving references to certain lattice 
data and to non perturbativc statements like solutions of SD equations 
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while many possibilities are opened for the renormalisation condition of the 
vertex. We need not specify it for reasons explained, below eq. (fTC|) . Then: 

9r{p) = gBGB{lJi^Y'^FB{tx^)IMp) (9) 

This imphes that the product gii{ii)zi{fi){G ii{k^ , fi))^^"^ Fii{k'^ , ^) is independent 
of /X. 

Now, gnip) would tend to a finite limit for small fi if eq. [3] would hold, 
under the additional assumption that zi(/Lt) is finite for — > ^. 

-The input of lattice data 

Recently, lattice data have also entered the game and have contributed much 
to the discussion, by showing features quite contrary to this consensus. Our 
motivation here is to try to clarify the situation within the SD approach by 
exhibiting new numerical solutions restoring the agreement between the lattice 
data, the numerical study in the continuum, and the analytical considerations 

6 

We can test the relation ^ on the lattice, by computing G{k)F{k)-^ , which 
according to the above relation ^ would be expected to tend to a finite value 
at small p. In fact, this is clearly contradicted by the lattice data, see Fig. 1 in 
our paper [S], and, from Sternbeck et al. Fig. 4 in [13], Fig. 3 in [IS], which show 
that the product decreases rapidly at small p, possibly to zero. In addition, 
as we will see from sections [2] and |3l for solutions satisfying the relation ([Hlthe 
product NcgR{fiyzi{fj,)GR{k^,iJ,)F]i{k'^,^)'^ must tend, when fc ^ 0, to a finite 
value which is much larger than the value observed at the smallest accessible 
momenta (it should be 10 tt^ in the case ac = !)• 

In addition, in the above studies, which try to solve the coupled SD equa- 
tions, the gluon propagator is also predicted, and it is found that ac is positive, 
which is anyway also suggested by lattice QCD. If so, and if the relation ^ 
would hold, then it would imply finally that ap < 0, i.e. the ghost dressing 
function should be singular. One could consider this result on the ghost dress- 
ing function to be theoretically appealing, because it meets another familiar 
statement, the Kugo-Ojima criterion for confinement (according to which the 
ghost dressing function should be singular). 

But, from the admitted values of ac > I, ap < —1/2, F(fc^) should present a 
power behavior close to 1/fc, or more singular. This stronger result is excluded 
by the lattice data, which allow at most a very weak singularity ; indeed 
according to our first analysis [5] and to the study of Sternbeck et al. [9] (see 
also their recent large volume study [lOj). the power seems at most ap = —0.2 
down to momenta around 0.3 GeV ^; in fact, we have obtained better fits of our 
own SU{2) and SU{3) data with logs rather than with powers [TT]; finally, if 
we abandon any prejudice, it appears that it is compatible with a finite value 
as well. Certainly, something close to = —0.5 is not possible, unless there is 

^Note that the UV finiteness of zi{fi) does not imply that is fi independent, even in 
perturbation, and in particular for the symmetric MOM scheme, see eq. 13 in our ref. [5], 
extracted from the results of Chetyrkin and Retey [^; however, we can suppose that the non 
perturbative IR behavior is not too singular. 

^An attempt to describe the lattice data within SD coupled equations is made in 0. For 
a recent attempt to accommodate the lattice data (with a finite non zero ghost dressing 
function) within the Gribov-Zwanziger approach, see [8] 

'^Note that these authors plot along the x axis in the Figures 
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a sudden change of behavior very near fc = 0. This conclusion is reinforced by 
the recent results of Cucchieri and Mendes at very large volumes [T2]. 

-Analytical setting 

In view of this situation, in our paper [5], we started a new discussion on 
the implications of the ghost SD equation for the IR behavior of the ghost 
propagator. 

We consider this equation in its subtracted, UV convergent form: 



GB{{q-kr)HMq.k) GB{{q-k'f)HMq,k') ] , 

k' is an arbitrary subtraction point, taken for simplicity parallel to fc, fc' = 

k\J fc'^/fc2. Hi is one of the invariants in the Lorentz decomposition of the 
ghost-gluon vertex: 

"^Blli-Q^k^q- k) = igBr'"'quTButi{-q,k:q~ k) 

= ^9Br"'{q^HiB{q,k) + iq-k)^H2B{q,k)) (11) 

where —q, k and g — fc are respectively the entering momenta of outgoing ghost, 
ingoing one and gluon, . 

In all our present considerations this equation is considered for given gluon 
propagator and vertex, and the ghost dressing function appears then as the 
solution to the equation. This is what we call the SD ghost equation. We do 
not try to solve any other SD equation. The advantage of concentrating on 
this equation is that it is much simpler than the gluon one or any other, to the 
point that analytical statements can be formulated for the ghost, with given 
gluon propagator and vertex. On the other hand, various assumptions on the 
IR behavior of the gluon propagator and the vertex may be used, in particular 
those advocated in the above references. Both these inputs and the output ghost 
solution can be tested through the lattice data. 

The renormalised form of this equation is : 



1 



GR{{q^kf)HiR{q,k) 
{{q-kW 



-Nc9rZi 



{k.qf 
k^q"^ 



GR{{q-k'f )HiR{q, k') 
{{q-k'W 



FR{q^) 



(12) 



where zi is the renormalisation constant of the ghost-gluon vertex. We know 
that, in Landau gauge, Hisiq, 0) + H2B{q, 0) = 1 which implies that zi is finite 
for any momentum configuration Let us remark that this implies that the 
subtracted SD equation is convergent. Indeed, zi and the l.h.s. of ea. p^ being 
finite the integral in the r.h.s. must be convergent. This was not obvious in the 
bare version. 



®Let us recall that what has been really demonstrated in the paper of Taylor |13j is the 
equation we have just written, i.e. for a vertex with zero ingoing ghost momentum. Then 
£l = 1 for this particular MOM renormalisation ; in general it will remain finite but different 
from 1. The detailed explanations on the Taylor paper are given in our article [5] 
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In the following, we set k'^ = to get the one variable renormalised integral 
equation. 

One can wonder whether the solutions of this subtracted SD equation are 
also solutions of the unsubtracted one : 



1 



(2^ 



1 



GR{{q-kf)H^R{q,k) 



This is seen to hold simply by making : 



FR{q^ 



(13) 



d^q 



1 + NcgiiZi 

{k.qf 



1 - 



GR{{q-kf) 



{{q-kW 



HiR{q,k) FR{q^) 



fc2 = 



(14) 



Of course, one has now to regularise the integral in some way, and this in- 
troduces a finite arbitrariness in Z3. The divergence of Z3, coming from the 
integration, will be of course cancelled, as it was in the subtracted form (fT2|) . 
by the divergence of the remaining part of the integral. 

-Introducing the regular solutions for the ghost dressing function 

We concluded in ref. (^) that, in general, under the usual IR regularity 
assumption for the ghost-gluon vertex, the SD ghost equation implies the rela- 
tion ^ by itself, without recoursing to the gluon equation. However, there were 
exceptions (see below), but we have first discarded them. Therefore, since the 
relation is definitely seen to be violated on the lattice, while the SD equation 
is automatically satisfied, we first suggested in the same paper, as a way out of 
this puzzle, that the vertex invariant Hi could be IR singular, instead of being 
constant. 

Then, it soon appeared, in view of the lattice data, in particular thanks to 
Sternbeck et al. - [IS] as well as to the previous work of Cucchieri et al. [IB] , 
that this possibility is very unprobable : indeed they measure HiB{q,q) (gluon 
at zero momentum, and contraction with g^), and they find it roughly constant 
and close to 1. Therefore, our attention has been drawn to the cases predicted 
in our analytical discussion of the SD ghost equation [S], where the relation [3] 
can be violated in spite of having a regular ghost-gluon vertex. These are the 
cases where ac — 1 and ap — 0, i.e. where the ghost dressing function is 
regular ^ at origin. As we have said, we did not pay attention to them in the 
beginning. But we have become aware that this possibility is attractive because 

1) on the lattice, aa seems not far from 1,[I3[TS|, i.e. the gluon propagator 
is not far from being finite (see also the very recent very large volumes stud- 
ies of the above references [l0],[12]). Thus it automatically leads to a rapidly 
decreasing G{k^)F{k'^)'^ , 0{k'^), behaviour when k approaches in agreement 
with what is observed. 

2) last but not least, on the lattice, the effective ap is compatible with 0, as 
we have seen above. 



^For some qualification of the term "regular" used in the present context, see the Intro- 
duction, below eq. 11511 
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Therefore, the appeahng possibihty ap — 0, has been adopted in our subse- 
quent paper [TT] , 

On the other hand, no statement can be made from these analytical consid- 
erations on the ghost SD equation as to which solution for the ghost propagator 
should be effectively preferred in real QCD. A complementary theoretical 
input comes from the Slavnov- Taylor (ST) identity for the three- 
gluon vertex. From this identity, we have demonstrated in that the ghost 
dressing function should be IR finite. 

The aim of the present paper is to reconsider this question by a numeri- 
cal study of the ghost SD equation with input from lattice data for the gluon 
propagator, and the simple and widely admitted constancy assumption for the 
vertex. The conclusion is striking : it is found that the IR finite solutions vi- 
olating the relation ^ indeed exist. Which type of solution is realized, either 
singular or regular at fc — > 0, depends on the actual value of the QCD coupling 
constant. One singular solution and only one is obtained, for only one value 
of which we call " critical" , and it can be completely exhibited ; according to 
our calculation, it cannot be the one of real QCD, because it disagrees 
grossly with the lattice results over a large range of momenta. There- 
fore, in agreement with our ST statement, the actual ghost dressing function 
must be regular, i.e. IR finite, and indeed we find solutions regular at fc ^ 
describing very well the lattice ghost data, with values of the coupling 
constant close to the one estimated from the actual bare coupling constant of 
the lattice. In summary, the combination of the numerical resolution and of 
the ghost lattice data or the ST theoretical input allows to discard the singular 
solution. 

-Warnings 

A caveat must be made now. To make the discussion simple and to keep it 
close to the commonly accepted conceptual framework, we have adopted above 
the usual assumption of a pure power IR behaviour for dressing functions. But 
this is by no means a necessary assumption; especially for a massless theory, 
it would not be unexpected to have a behavior with log factors accompanying 
integer powers of fc^ . For instance, according to our Slavnov- Taylor discussion 
[SI 1191 I20j . the gluon propagator must be infinite at the origin, under some 
regularity assumption for the three-gluon vertex. Then a way to reconcile this 
statement with the observation of an apparently IR finite propagator from the 
lattice is to assume that this infinity is logarithmic, which would make it very 
difficult to detect on the lattice The behavior of G{k) would be G{k) ~ 
fc^(log(fc^))'' {h' > 0). We present the results disregarding the logs. We have 
also checked in our numerical calculation that including such a log in the 
gluon propagator does not change appreciably the ghost propagator 
deduced from the SD equation. 

In addition, in our analytical discussion [19] , we have shown that if ac — 1 
and ap = 0, one must have in the ghost dressing function logarithms of the 
type 

F{k) =a + b k^\og(k^). (15) 

"'^''Of course, the gluon propagator is always finite at fc = in a finite volume. What we 
mean here is that the gluon propagator at fc = seems to be relatively constant for increasing 
volume. And, one cannot exclude a very slow variation, such as expected from a logarithm. 
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The effect of such logarithms is very weak, so that for the present purpose we 
count this as "regular" , and anyway, it is IR finite. Nevertheless, we are able to 
display the effect of this logarithm in the numerical calculation. 

As we have said, we have checked that including a log in the gluon propa- 
gator does not change appreciably the ghost propagator deduced from the SD 
equation. This is in agreement with an analytical discussion, which shows that 
only the power of the log in eq.dTHl) is changed. 

-Assumptions and inputs of the calculation 

To summarise, our starting assumptions for the numerical calculation below 
are finally the following : 

1) we take the ghost gluon vertex invariant Hi{k.q) implied in the SD ghost 
equation (see below) as momentum independent. This is a rather usual assump- 
tion made in SD studies, and as we said above, it is in rough compatibility with 
present lattice data. In fact, it would be sufficient to assume simply a regular 
vertex to get the same qualitative conclusions, but for a numerical study we 
have to make a definite choice. Anyway, let us emphasize that our first goal is 
not to make a realistic quantitative prediction, but to demonstrate that, even 
with this type of regularity assumption, one can obtain solutions not obeying 
ug + 2aF = 0, contrarily to the common belief. For this purpose, it is not 
required to bother about what would be the most realistic assumption for the 
vertex. 

2) for the gluon propagator, in the small momentum region, we use an in- 
terpolation of the gluon propagator given by the lattice, with ac — 1. 

We are aware that several SD studies, for example [I], exclude this latter 
possibility when considering the coupled equations, since ac and ap are then 
determined separately, and they found then that k = —ap = l/2aG > 1/2, 
therefore ac > 1. At this point, we recall that the paper of Bloch [3] concerning 
the coupled equations finds solutions with ac = 1 and ac + 2aF = {k = 
1/2), thanks to a more refined treatment of the gluon SD equation. He can 
then reproduce rather well the gluon lattice propagator. Then the question is 
compelling : knowing that the gluon lattice data are well reproduced by him, 
and that the lattice data satisfy the ghost SD equation he is solving, how can it 
be that the lattice ghost dressing function differ from his prediction 1/fc ? We 
show that this is not due to lattice artefacts, but to the neglected possibility 
that there are different type of solutions to the ghost equation for the same 
gluon propagator, depending on the value of the coupling constant. 

2 Analytical considerations on the behaviour at 
small k 

As a preliminary to the numerical study, let us recall or establish analytical 
relations, which can be used as tests of the soundness and accuracy of our 
numerical calculation. We extract them from a more complete discussion, which 
will be given in a forthcoming paper 21]. 

In this section and hereafter, since we adopt the constancy assumption for 
the ghost-gluon vertex, it appears immediately that the coupling constant only 
appears in the combination : 

EE N,g%ziHiR - N^glZ^Zj/ z^Hir = N^glZ^Z^HiB (16) 
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where Hij^ or Hib are constants. We use this auxiUary notation througout the 
rest of the article. From the last equality, it is obvious that is independent 
of the way one renormalises the vertex, since only the bare vertex takes place. 

We proceed as in our paper [S] : - we separate the integral into one UV part 
and an IR part, respectively for q > qo and g < go- 

-for the infrared contribution (see eq.(14) and (17) of [3] with ar = and Hi 
and h constant), defining A and B by FR(fc) ~ Aik"^)"" and GR{k) ~ B{k'^Y° 
when /c ^ 0, i.e. for k < qo. 

-we write eq. (|12p replacing k by Afc, taking k' — Xnk and performing the 
change of variable q Xq. We then consider the IR limit A ^ 0. 

1) Singular solution. Let us establish a relation for g^, the value of 
corrresponding to the singular solution, which is found to be unique. It is 
inspired by Bloch [3]. 

One can write, at leading order in A : 



(A2fc^)--(i - --'-n - -5^(A^)-+"-^^i? 1'^*^' ('^')"'^ (i - 

X [{(q-k)'r^-^-{{q-Kkfr--'] 

This integral being C)(A^^"^+""^^) dominates over the UV part which is 
C'(A^), by a negative power of A, for aa = 1 and ap < 0. We can then 
neglect the UV part. This gives the relation 2aF + ctG — 0, whence ap = —1/2, 
Fiji^) ~ 1/fc. Moreover the integral of the r.h.s. can be analytically performed. 
Defining the function: 

1 V{2 + a)T{2 + h)T{-a-h-2) 
^) - 16^2 r(-a)r(-5)r(4 + a + 5) ^^^^ 

its value is equal to (1 - K"^"^)(fc^)~"^<i>(aG) where: 

$(aG) = 

-i(/(-^,«G - 2) + /(-^,aG - 1) + /(-^ - l,aG - 1)) + 

- 1, - 2) + f(-^ - 1, cg) + /(-^ + 1, aG - 2)) (19) 

leading to 'g^A^B = In our case, q;g = 1 and ap = -5, $(1) = 

and the relation becomes: 

52a2G^2)(0) = 10^2 (20) 

(2) 

where G]^ is the gluon propagator. 

This is only a relation between 'gl and A and it doesn't allow us to know 
a priori the value of 5^ before any numerical computation, unless a very small 
renormalization point ji is choosen. In this case we have: A = (/x^)~"^ and 
B = (^2)-"G so that gl = ^ {gl = 10^^ ^^^^ ^ 

2) Regular solutions. If, on the other hand, ap — then the l.h.s. is 

trivially zero at leading order in A, and one has to go a step further in the 

(1) 

expansion to get a non trivial result. Noting that A = Fji{0) and B = G]^ (0), 



8 



where G^^^ is the gluon propagator, finite at the origin, the IR part of the 
integral has the form : 

-ff'^A2fc'(l - log(go/(A fc))F«(0)2Gg)(0) + ... (21) 

where the dots denote subleading C'(fc^) terms. It still dominates over the UV 
part, although this time it is only by a logarithm. We then write consistently 
Fr{P) = a + 6fc2log(l/fc2) + 0{P) in the l.h.s, and we find: 

Fnie) = Fn{0) (^l -g^^Fn{OrG^^\o)eiog{Mye)^ (22) 

being some scale which we cannot derive from this IR expansion. 



1 1 

{q — kY (q — nk)^ 



3 Numerical solution of the Schwinger-Dyson 
equation for the ghost 

In this section we want to see if the two types of solutions {ap — and 2aF + 
OLQ — 0) suggested by our analytical discussion in [3 [H] actually exist for the 
same gluon propagator. We answer positively by solving numerically the ghost 
SD equation for given gluon propagator and vertex. In the following we shall 
use the subtracted form of the Schwinger-Dyson equation for the ghost in the 
Landau (i.e. Lorentz) gauge. This equation has been written above, eq. ()12p . 

We start from an IR finite gluon propagator (ac = 1) extracted from our 
lattice data in pure Yang-Mills theory, with Wilson gauge action, /3 = 5.8 and 
a lattice volume equal to 32'*), for momenta lower than 1.5 GeV ; this choice is 
justified to have moderate UV artefacts. We extend it to larger momenta using 
a one loop asymptotic expansion (with h^MOM = 1 GeV corresponding to the 
standard ^j^g = -240 GeV of lattice quenched QGD). On the other hand, we 
take Hi{q, k) to be constant with respect to both momenta As we said above, 
this is suggested by the lattice data for q — k (i.e. for zero gluon momentum), 
but we extend it to all values of q and k. The authors of ref. T5] find a bare 
vertex very close to 1 in this zero momentum gluon configuration for a large 
range of 

We work in the MOM scheme, and set fc'^ appearing in eg. p^ as the 
squared renormalisation scale /^^ (/i has been chosen at an optimum 1.5 GeV, 
not too high to allow the lattice data to be safe, and not too small to display 
the differences between solutions at small momenta). The equation we have to 
solve is : 

1 



i^fi(F) " J (27r)4 V 

GR{{q-kf) GR{{q~k'f) 
{{q-kfY {{q 



(23) 



^^This cannot be an exact statement, as already shown in perturbation by the calculations 
of ref. |22l |6] ; although finite, the vertex invariants do depend on the momenta through the 
running as 
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Note that this equation implies that depends only on the renormalisation 
point chosen for the propagators ; it is independent of the particular way used 
to define the renormalisation of the vertex, in agreement with ea. (|16p . Eq. [^3] 
can be still transformed to a new form which makes the numerical calculation 
and the presentation of the various solutions easier ; for this, we subtract the 
equation at fc = 0, to let the value of Ffi{k) at the origin appear and to eliminate 
the reference to the particular normalisation point /it, and we redefine also the 
unknown function to be calculated as F{k) = gFif{k). Then the reference to 
the value of 5 also disappears; we end with : 



1 



1 



F{k^) 



m 

GR{{q- 
{{q^ky 



d^q 



1 - 



{k.qf 



(27r)4 V 



(24) 



We solve this equation ((^ for F{k'^), for a set of values of F(0). It is easy 
to see that from this solution we can reconstruct the desired solution of eq. ((23|) 
for any renormalisation point and any value of Indeed, we have ^(/i) = F{^'^), 
so that for given /i and 5, we have just to identify the value of F{Q) such that 
^^(^2) ^ ~ -Y\ie\i, we reconstruct FR{k'^) through F^Xfe^) = F{k'^)/g{fi) 

By construction, all the solutions found in this way are finite at the origin. 
The solution divergent at the origin will be found by setting -^j-^ = in eq. 

([24|) . It can also be approached by making F{Q) larger and larger. 

We have looked for solutions of eq. (f24|) with the integral cut in the UV at 
q ^ 30 GeV. We have discretized it in k and q. Taking values of the momenta 
spaced out by 0.01 GeV for g < 2 GeV and by 0.1 GeV for q > 2 GeV we have 
computed the angular integral of the r.h.s. of eq. ((^5|) . Then, we solved this 
equation by iteration. Minus the integral in the r.h.s. is positive, allowing an 
easy convergence. We linearize it at each step, following the Newton method, 
to accelerate the convergence of the iteration procedure, as suggested by Bloch. 

The results are the following: 

1) Critical case, singular solution. We find a solution with 
i.e. F{0) = 00. We find then the corresponding "critical" constant: 

g^^ = ^^(1.5 GeV) = 33.198.... 

The relation of eq. ([20|l happens to be very well satisfied: 

g;'^'G(2nO)T7^ « 0.994.... 



= 

F(0) 



(25) 



(26) 



The integration near k = can be improved by taking explicitly into account 
the analytical behavior of the kernel, and assuming that the solution behaves as 
1/k at smah k. This imposes eq. and one checks that g^k'^F{k'^)'^G'-^\k'^)j^ 
goes very smoothly to 1 when fc ^ 0. 

2) Regular case. We find a solution for all F{0) > 0, and only one for 
each F{0) > with our method of solution. . From our numerical solution at 
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^ ~ 29, which corresponds to the best description of lattice data (see Fig. (2), 
we can test the relation giving the fc^ log(fc^) term. The result is presented 
in FiglU and the slope agrees well with what is expected from ([2^ : 4.06 against 
4.11. 

The critical value of the coupling constant, as well as the corresponding curve 
of F(k), can be very well approximated by the regular solutions at very large 
F{Q). When F{Q) is larger and larger, the eq. is valid only in a smaller 

and smaller region near fc = 0, while in an intermediate region, we observe a 
1/k behaviour. 




Figure 1: The a + bk^ log(fe^) fit at small momentum (dashed line) to our continuum 
SD prediction for the ghost dressing function, renormalised at /x = 1.5 GeV for = 
29. (solid line) ; the slope of the fe^ log(fc^) term is 4.06 ; the agreement with the 
expected coefficient of log(fc^), 4.11 from the eq. (|22|l . is striking 



In conclusion, in the case aa = 1 we have exhibited a continuum set of 
IR finite solutions for arbitrary -F(O), and a unique singular solution for = 
g^, with ap = These are the only solutions obtained with our iteration 

procedure using the Newton method. Of course this doesn't prove that no other 
solution exists. 

3.1 Why the regular solutions have not been obtained pre- 
viously 

The regular solutions could not be obtained by the proponents of the equation 
131 because, as it seems to us, they discard them from the beginning, and thereby 
choose the critical value of the coupling constant, by making an implicit assump- 
tion when solving the so-called "infrared equation" for the ghost SD equation. 
One can see this in the papers by R. Alkofer et al. (for instance [I], eqs. (43), 
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(44)), or in the detailed discussion of Bloch |3|), eqs. (55) to (58) . 

Let us explain this briefly. They consider the above unsubtracted equation 
(note that this requires then an UV cutoff, which we avoid in our discussions 
by considering the subtracted form-see below, next section-) ; we write again 
the unsubtracted form : 

- - NcOrZi / — — 1 



Ffl.(fc2) "^^-^ 'J (27r)4 V fc2q2 

-GR{{q^k)^)Hin{q,k) 



{{q-kYf 



FR{q') (27) 



One tries to match the small k^ behaviour on both sides of eq. [27] This is 
done for example in eq. (58) of [S]. A condition is then written which consists 
in equating the coefficient of (fc^)""'' with the corresponding one in the r.h.s.. 
However, one notices that on the r.h.s.., there is a constant contribution a {k"^)^ . 
Therefore unless the constant term Z3 is cancelled by the integral contribution 
for A: — > 0, we have necessarily ap = 0. To have ap < as the author finds, one 
needs this cancellation. This is what is implicitly assumed, but not stated 
explicitly. The condition of cancellation is : 



(28) 

fe=0 



However, this additional equation does not derive from the starting SD ghost 
equation, and indeed it is not satisfied in general by the solutions of this basic 
equation, as we show by displaying actually IR finite solutions. In fact, it can 
be valid only for a particular value of the coupling constant, the critical 
one which is solution to the equation of Bloch, his eq. (58), and which we derive 
rigorously through the subtracted equation (see our eq |20p . 

Let us make precise that at this stage the value of the coupling constant is 
taken as a free parameter. It should be however fixed in the end consistently 
with the lattice data which we are using, which is done in the next section. 



4 Phenomenology 

Having presented the general study and found the announced two types of so- 
lutions for the ghost dressing function, either regular or singular, the question 
is : which one is effectively realised on the lattice, and therefore in true QCD ? 
Let us recall the classical problems which hamper the answer : it is not possible 
to know from the lattice data with total certitude whether the ghost dressing 
function is singular or not, because 1) on the one hand, the "singular" qualifi- 
cation in itself does not tell how close one should be to the zero momentum for 
the singularity to show up ; 2) on the other hand, one cannot get arbitrarily 
close to zero momentum on the lattice. 

A better mean is offered by our calculation : it predicts the behavior of the 
respective solutions for the ghost over the whole range of momentum, and 
not only very close to fc = ; then looking to the lattice data, we can identify 
which is the more compatible with the data From fig. we see that we can 

-"^■^At this stage, it is useful to stress the advantage of working with the renormaUsed form 
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discard rather safely the singular solution, because it is passing much above the 
lattice data points over a very large range of momentum: around 50% at 
the leftmost point measured on our lattice, k = 0.26 GeV, but still quite sizable 
near k — 0.5 GeV. It is therefore quite unprobable that any lattice IR artefact 
could fill the gap. The advantage of our method is that, by calculating what 
the critical solution should be at rather large momenta, we are able to discard 
it more convincingly. 

On the other hand, we find a very good description of the lattice data in 
the range = 28.3 — 29.8 (the range is defined by one standard deviation 
except for the lowest point). This striking agreement is illustrated by Fig. [2] 
(For indication, we quote the IR limit Fii{0) = 2.51 for the same fj, = 1.5 and 
= 29). 



3 \ 

\ X dressing function lattice result 

_ \ best SD solution 

\ singular SD solution 

2.5 - \ 




I I \ I I I I I \ I I 

0.25 0.5 0.75 q 1 1.25 1.5 

Figure 2: Comparison between the lattice SU(3) data at /3 = 5.8 and with a volume 

32* for the ghost dressing function and our continuum SD prediction renormalised 

at fi — 1.5 GeV for ff^ = 29. (solid line) ; the agreement is striking ; also shown is 

the singular solution which exists only at = 33.198.... (broken line), and which is 

obviously excluded. 

Moreover, we can perform the following consistency test. We start from the 
equation (|16p defining g^ in terms of bare quantities and we apply it to connect 
the continuum g^ to the lattice bare quantities: 

g2 = N.glli = N,{6/p) Fl,{ii^)GB{f^^)HiB (29) 

of the SD equations ; indeed the continuum and lattice versions are more directly comparable 
than the bare ones. As we have seen in our paper 0, the bare lattice equation for the ghost 
is affected by an important artefact which vanishes only very slowly with the cutoff, being 
of order O(g^). In the renormalised version, this effect is included in the renormalisation 
constant Z3, and we are left only with the much smaller cutoff effects of the type 0{a") 
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We then ask whether our range = 28.3 — 29.8 is reasonably consistent with 
the r.h.s. of ea. iP^ as given by lattice data. Let us stress that eq. ((^^ should 
be then only approximate in several respects; first, it is valid up to finite cutoff 
effects, as well as volume effects ; second, we have replaced the lattice vertex 
invariant HiB{q,k) by the constant His, which is very rough ; and we cannot 
test on the lattice the assumed constancy over the momenta which are actually 
implied in our calculation. We have only at our disposal a lattice measure of 
HiB at q = k. Last but not least, lattice measurement of such vertex quantities 
is difficult; it is very noisy. Then, the test is only qualitative; nevertheless, the 
result is very encouraging, as we see now. 

Indeed, from the above value of found in the continuum on the one hand 
and, on the other hand, from the lattice data /3 = 5.8, G'B(/i^) ~ 2.89 and 
~ 1.64 {fi is here chosen as 1.5 GeV), we can deduce the value of the 
factor HiB needed to satisfy equation (j29p . and which represents some average 
on momenta. We find Hib — 1.2. This number should be compared to the 
lattice measurements which are for Hib at g = fc, and which give about l.,with 
large errors (slightly larger or equal to 1., according to the k value, see ref. [15j). 
The comparison seems very encouraging considering the large uncertainties of 
the procedure : lattice artefacts, errors on Hib, and finally the fact that our 
HiB is some average on momenta away from q = k. 

Another way of presenting the striking difference between the regular solu- 
tion and the singular one is in term of the familiar product discussed in the 
introduction : G nik) Fii{k)^ From the analytical discussion, in the critical 
case, it should tend to 10 jg^ — 3. when /c — *■ , and numerically it should 
be 3.14 at our smallest lattice momentum k = 0.26 . This is completely at odds 
with the lattice : the lattice value is 1.28 at k = 0.26, with a clear tendency to 
still lower values at smaller k. On the contrary, our regular solution fits per- 
fectly the lattice data. We illustrate this in Fig. [31 Let us stress that our SD 
solutions (continuous curve) are obtained in the continuum and in infinite 
volume ; they appeal to the lattice data only to have a physically reasonable 
definite gluon propagator as input to the SD equation. 



5 Conclusion 

The relation ac + 'iap = is usually believed to be an unavoidable consequence 
of SD equations. The problem is that the lattice data grossly contradict this. 
Indeed, for small momenta G(fc^)i^(fc^)^ — *■ very fast. We resolve this con- 
tradiction in the following way. We show that this belief is wrong and that an 
alternative exists by solving numerically the ghost SD equation with input from 
the lattice for the gluon propagator and the vertex . 

The alternative is the one we have previously envisaged as a possibility in a 
general analytical analysis OHI]: ar = withac ^ 1, therefore G(/c^)F(fe^)^ ^ 
as shown by the data, but in the present article, the existence of such a solu- 
tion is demonstrated by actually solving in F the equation for a given G. This 
solution violates the statement ac + 2aF — 0- 

Usually, this quantity is presented with multiplication by an additional factor including 
the renormalised coupling constant and possibly other factors. Here, we present the raw 
product to avoid any ambiguity in such procedures. 
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0.0001 0,001 0,01 0,1 1 10 

Figure 3: Comparison between our lattice SU(3) data at /3 = 5.8 and with a vol- 
ume 32* for the product of gluon times ghost square dressing functions GR{k)Fii{k)^ , 
renormalised at fi = 1.5 GeV, and the corresponding curve for the continuum singular 
solution ao + 2qf = 0, which exists only at = 33.198....; it is obviously excluded; 
up to a factor /A-k, it corresponds to the a(fc^) presented in Fig. 8 by Fischer [23j : 
the form is very similar. Also shown is our continuum regular solution for — 29. 
(solid line) for which the agreement is striking. 



The relation uq + 2aF = would imply that the ghost factor is singular, 
since ac > 0. The numerical solution of the equation then adds another strong 
reason for rejecting this relation. A singular solution - which necessarily satisfies 
ac + 2aF = - only exists for a definite value of the coupling constant. We 
calculate it and find that it grossly differs from the lattice data on a large range 
of k, and not only for the smallest momenta. 

The alternative solution which is regular (a^? = 0), is realised if the coupling 
constant is smaller than a certain critical value, while the singular one is present 
only at the critical value. The lattice data for the ghost are very well reproduced 
for a coupling constant close to the one expected from the value /3 = 6.0 used 
for the lattice calculation. We are then confident to have found the actual QCD 
solution, up to moderate artefacts. 

Our numerical study therefore adds strong new arguments from the lattice 
data in favor of this alternative . It is based on the ghost equation only, since we 
feel that the gluon equation, being much more complicated, suffers from much 
more uncertainties, due to the necessary critical approximations to be made - 
this has been illustrated by the findings of Bloch [5] . 

The alternative, regular, solution has not been found in usual studies, be- 
cause they have chosen by construction the critical value. 

The important physical consequence is that we do not get the alleged non 
trivial IR fixed point for the MOM coupling constant since gnik) when 
fc ^ at fixed go as for the three-gluon couplings, in agreement with lattice 
data. At this point, it is important to insist on the fact that there are infinitely 
many definitions of "the" QCD coupling constant. A priori, there is no reason 
for their IR behaviour to be universal. In particular, the ones defined from 
elementary fields Green functions have no reason to behave in the same way 
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as more physical definitions such as taken from the pcrturbative expansion of 
certain physical hadronic amplitudes. 
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